Green’s theorem 1 


Chapter 12 Green’s theorem 


We are now going to begin at last to connect differentiation and integration in multivariable 
calculus. In addition to all our standard integration techniques, such as Fubini’s theorem and 
the Jacobian formula for changing variables, we now add the fundamental theorem of calculus 
to the scene. In fact, Green’s theorem may very well be regarded as a direct application of 
this fundamental theorem. 


A. The basic theorem of Green 


Consider the following type of region R contained in R?, which we regard as the x —y plane. 
We use the standard orientation, so that a 90° counterclockwise rotation moves the positive 
x-axis to the positive y-axis. Then we assume the existence of two continuous functions a(y) 
and b(y), defined for c < y < d, where a(y) < b(y) for c < y < d and 


R= {(x,y) | aly) < x < bly)}. 





x 

















Given a function R —> R of class C! , we then compute by means of Fubini’s theorem 


b(y) 
/ a = [ [. or dey 
Ox 
R 
PrC r=b(y) 
2 [ F(2,y) 


dy 
_ [ Few.nay— f F(a(y), y)dy- 


v=a(y) 
We now write the right side of this equation as 


/ Fdy, 
bdR 
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where we mean by this notation the counterclockwise integration of F’, restricted to bd R, with 
respect to y. (The actual definition is given in the formula.) 


bi 


Notice that on a horizontal portion of bdR, y is constant and we thus interpret dy = 0 there. 


B. Line integrals 


We have now met an entirely new kind of integral, the integral along the counterclockwise 
bdR seen above. Before proceeding further, we need to discuss this sort of oriented integral. 

It will prove useful to do this in more generality, so we consider a curve y in R” which is 
of class C!. Thus, [a, b] —, R” is of class C!. We typically denote the independent variable 
(the “parameter” ) as t. If f is a continuous real-valued-function defined on R”, we then define 
the line integral 





























[see = [ towne 


Here of course 1 <i <n. The notation is intended to be very suggestive and to lead us to 
the appropriate substitutions z = y(t), x; = 7:(t), and dx; = 7;(t)dt. 

It is an easy matter to imagine some useful properties of this sort of integral, and even 
easier to prove them. 


INDEPENDENCE OF PARAMETER CHANGE 

Suppose first that y(t), a < t < b, is represented instead as y(h(s)), c< s < d, where h is 
a C! function such that h(c) = a and h(d) = b. Then using the parametrization y(h(s)) leads 
to the line integral 


J toms Zrtisas =f) Fevln(s)))aflats))H(s)as 
OF powmyiwa, 
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which is the original line integral. ‘ 





This proves the desired independence. 
On the other hand, if instead h(c) = b and h(d) = a, then we obtain 


b 


[ tomy Zatisas = = [ro mynitoee 


a 


so we get the anticipated change of sign. t 





LINEARITY 
This is virtually obvious from the definition: 


[ofan = af fae, if a is a constant; 
¥ 


[r+ oe — [ties f 9 


COMBINING CURVES 

Suppose we are given curves y and 6 such that the final point of 7 equals the initial point 
of 6. Then we can think of a new curve that covers y and then 6. This curve will be continuous 
but not necessarily C!. We denote this curve loosely as y + 6, and then the clear result is 


: je; = [tans f faz. 
+6 7 6 
GENERALIZATION 


If fi,.-.., fn are given functions, then we define 


J fides +++ fader = [titer tent f fade 
oy Y Y 
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Another way to think of this is to arrange the functions into a vector F = (f,,...,f,) and 
formally write 


> 


ES (eis cigalig 
dx (Wii, ca2g hy), 
Fed&t = fidx,+---+ fnrdrty. 


Thus we have the abbreviation 


[Pea 
y 


FUNDAMENTAL THEOREM OF CALCULUS 
Suppose that f is a real-valued function of class C' and that 


F = Vf = (f/dm,...,0f /Otn). 


[fae +---+ f,dtn. 
- 


Then we calculate directly 
b 
[vie = i, Vi(y(t)) ey (t)dt 
7 a 
er, 


chain rule 
mene fF l(t) 





Thus we have proved 
[ve edz = f(final point of y) — f(initial point of y). 
- 
We shall say much more about this equation in Section E. 


C. A general Green’s theorem 


We now return to the formula of Section A, 


/ | or dey = / Fdy. (x) 
R 


bdR 


Green’s theorem 5 


The right side is now completely understood as a line integral taken along the curve bdR with 
its counterclockwise orientation. We presently have severe restrictions on what the region R 
can be, and we now show how we may easily remove many of these restrictions. 

Rather than give precise definitions to delineate the allowable regions, we prefer to rely on 
pictures. Suppose then we have a region in R? that looks like this: 














¥ 


$f, 


Such a region fails to satisfy the conditions of Section A. However, the desired equation (*) 
still holds for it. An easy way to see this is to divide R into two regions R, and Ry by means 
of a “cut,” in such a way that (*) is valid for each piece separately. 


¥ 


$b, 


We thus have 


Fdy, 


OF 
—dad 

// an nw 

Ry bdRy 

[f Spaeet = / Fdy. 
Ox 

Ro 


bdR2 


Now add these two equations. The left sides produce the double integral over FR itself, and the 
right sides combine in an interesting way. Namely, the part of each line integral corresponding 
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to the cut appears in each, but with opposite orientations. Thus they cancel to produce just 


[Fw 


bdR 


We therefore find that («) remains valid for this R. 


Of course, this procedure may be extended to several cuts, so that (*) remains valid for a 
region such as 





One more generalization allows holes to appear in R, as for example 





We again make cuts to get a couple of regions 
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Then we apply (*) to Ry, and Rz and add the results, noting the cancellation of the integrations 
taken along the cuts. The result still is (x), but with an interesting distinction: the line integral 
along the inner portion of bdR actually goes in the clockwise direction. 

A convenient way of expressing this result is to say that (*) holds, where the orientation 
of bdR is such that the interior of R is located to the left of the path of integration. 

Finally, we close this section with the corresponding result 


[ss —dady = - f Pas. (x) 


bdR 


The minus sign appears naturally, as we see from rewriting the basic case of Section A. Another 
way to think about it is to realize that the coordinate system with y first, x second, has the 
opposite orientation to the one we have been working with. 

Most texts combine these two formulas into a single one by using different letters for the 
two cases, and adding. Thus, 


|| (2-5) dady = f Par+ Qay 


bdR 














D. Areas by means of Green 


An astonishing use of Green’s theorem is to calculate some rather interesting areas. These 
come for example by using F' = x in (*) in Section C to obtain 


area(R) = [vv 
bdR 
Likewise, using F = y in (**) produces 
area(R) = =f yee 
bdR 


The latter equation resembles the standard beginning calculus formula for area under a graph: 


M4 
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[ feo 
= f ya 


It is interesting that the sum of these two formulas is often more easily exploited. Or use 
Q =x, P = —y in the formula at the end of Section C. The result is 


The integral 


is exactly the line integral 


taken around bdR. 


1 
area(R) = 5 / xdy — ydz. 
bdR 


This is often useful in situations where x and y appear rather symmetrically. 


ELLIPSE. We of course know from several investigations that the area inside the ellipse 
x? /a* + y?/b? = 1 is equal to tab. This example gives a nice illustration of our new formula. 
Namely, parametrize the ellipse by x = acost, y = bsint, where 0 <t < 27. Then 


xdy — ydx = acost- bcostdt — bsint - (—asintdt) 
= ab(cos’t + sin? t)dt 
= abdt. 


1 Qn 
= | abdt = zab. 
2 So 


Don’t miss the fact that if we had just integrated xdy or —ydz to find this area, the final 
integration would have been significantly more difficult. 


Thus the area equals 


FOLIUM OF DESCARTES. Here is a truly significant example, one for which our formula 
provides what is perhaps the easiest method for calculating a certain area. We first discussed 
this curve in Section 6E, where we used the implicit equation x? + y? = 3zy. In terms of the 
ratio t = y/x, we have the parametric presentation 


Y 
3t 
Ll = -5.—} 
14+ 28’ 
3t 
Yy = 


1+ 


Green’s theorem 9 


Now we find the area of the “leaf,” that region enclosed by the part of the curve in the first 
quadrant. Here the parameter t runs between 0 and oo. We calculate 


zdy —ydx = xd(tx) — tadxr 
x(tdx + xdt) — tadx 
= adi 

9t7 dt 
tera 


—3 
= a( 5). 
1 f* —3 
area = s/f i(sSs) 


a 1 


2148 
3 


Thus 








0 


PROBLEM 12-1. Obtain the area of the leaf of Descartes by using just the formula 


area(R) = = f yee 


PROBLEM 12-2. Show that the equations 


area(J) / ydx 


ady 
bdR 


can be regarded as equivalent by means of an integration by parts. 
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PROBLEM 12-3*. Find the area between the folium of Descartes and its asymptote 
x+y = -1, as shown in the figure. 


PROBLEM 12-4. Find the area enclosed by the curve 


in the first quadrant. 


PROBLEM 12-5. Consider the curve defined by the equation 2* = ry? + 4°. 


a. Use the parameter t defined by y = tx to express the curve in parametric form. 
b. Show that this curve has a “leaf” in the fourth quadrant. 
c. Calculate the area of the leaf. (answer: 1/210). 





E. Conservative vector fields 


Now we return to our discussion of line integrals in general, as they were introduced in 
Section B. 


DEFINITION. Let D bean open subset of R”. A vector field on D is a continuous function 


D +3 R’. 





Thus F assigns to each point « € D a point F(x) € R”. We always think of F(x) as a vector 
attached to the point 2. Thus we have a picture of the following sort: 
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F (x) 
* F (x) 


Using the components of F’ = (F\, Fb,..., Fi), we have defined line integrals of the form 


y 


In Section B we introduced the notation 


i, Fedzx 
af 
for these line integrals. 
Now we state an important basic theorem about these vector fields. 





THEOREM. Let D be an open subset of R", and let F' be a vector field defined on D. Then 
the following three conditions are equivalent. 











1. The line integral [| F edz is independent of path, in the sense that its value depends 
only on the initial point of y and the final point of y. 


2. The line integral a e dx = 0 for every loop in D (a loop is a curve y whose initial 
point equals its final point). 





3. There exists a function D .R of class C! such that 


F = Vf. 
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PROOF. The proof that 1 => 2 is just based on the simple observation that a loop which 
starts at xq and ends at x and the constant path which just stays at x9 are two paths with the 
same initial and final points. As the line integral over the constant path is zero, condition 1 
implies that the line integral around the loop is zero as well. 

The proof that 2 = 1 is based on another simple observation: if y, and y2 are two paths 
with the same initial and final points, then the path that goes along y, and then backwards 
along 72 is a loop. Thus condition 2 implies 


[Peas | Fedz = 0. 
1 y2reversed 


This proves condition 1. 

The proof that 3 = 1 is an immediate consequence of the last formula in Section B. 

The really significant part of this theorem is the proof that 1 => 3. 

Before we introduce the crucial construction of the function f we notice that we may as 
well assume that our open set D is connected. In other words, that any two points in D can 
be joined by some curve lying in D. For we may work on each connected component of D and 
produce the required function f on each of them in turn. 

In this context it is important to observe that any two points in the connected set D can 
in fact be joined by a “nice” curve, say one of class C! or one which is polygonal. For we need 
to be able to do line integrals along such curves. 


D 





Now our definition of f is actually forced upon us. For given the vector field F’, if there is 
to be a function such that F' = Vf, then we know from Section B that 


f(final point of 7) — f (initial point of y) = fF edz. 


7 


Now choose an arbitrary but fixed point b € D. It will serve as a “base point.” Then for any 
x € D we produce a curve y in D whose initial point is b and whose final point is x. Then 
our “potential” function f has to satisfy 
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Since b is fixed, f(b) is a constant; we may as well call it 0, as adding a constant to f does 











not disturb our desired equation Vf = F’. Thus we now define D oR by the equation 


fe) = Peas, 





where y is any “nice” curve in D starting at b and ending at x. 

We have now in fact used the hypothesis 1, since the value given for f(x) does not depend 
on the choice of y. 

Finally we must verify that Vf = F’. We’ll present two proofs. These proofs are indeed 
related, but they have distinct emphases. 


FIRST PROOF. Consider any fixed point y € D. Let y be a “nice” curve in D starting 
at b and ending at y. Also let B(y,r) be an open ball contained in D. For any unit vector 
h € R” and any 0 <s <r we use the curve 7 followed by the curve which is the line segment 
from y to y+ sh. Then 














fy +sh)= [ Pear | Fedi 


Yo [yytsh| 


= constant +f F(y +th) e hdt. 
0 
Yo 


B(y,r) 


We immediately obtain from the FTC the formula for the directional derivative 
d : ee 
qs fy + sh)) = F(y+ sh) eh. 
Set s = 0 to achieve . : 
Df(y;h) = Fly) @ h. 
In particular, when h = unit coordinate vector é;, 


Of 
Ox; 





(y) = Fi(y). 


Thus the partial derivatives Of /Ox; are equal to F; and are thus continuous; we conclude that 
f is of class C! and Vf = F. 
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SECOND PROOF. This time we go straight for the definition of differentiability as found 
in Section 2E. We have of course a candidate for V f(y), namely Fy). So we are led to 
consider the quantity 

fy +h) —fy)—Fly)eh 


for small h € R”. Then the definition of f gives 


Fly +h) —F0)—Fen= f hte Fe 





= [Fy +tnyenat— ry) en 


= / ‘FW + th) — F(y)]  hdt. 


Given € > 0, the continuity of F at y implies that there exists 6 > 0 such that for ||z —y|| < 6 
we have ||F'(z) — F(y)|| < «. Thus we have for ||h|| < 6 the inequalities 


f(y +h) — f(y) — FQ) ehl < / IF(y + th) — F(y)] © hlat 


< | |F(y + th) — F(y)|| |hllat 
Schwarz Jo 


1 
< i el|h|lat 
0 


= eh]. 


This is precisely what we needed! We conclude that f is differentiable at y and that V f(y) = 
Fy). QED 


DEFINITION. A vector field which satisfies the above equivalent conditions is said to be 
a conservative vector field. And any function f satisfying 3 is said to be a potential function 
for F’. 

That theorem is really wonderful: clean statement, easy but significant proof. However, 
it might be difficult to apply in practice. How would we ever verify the hypothesis 2, for 
instance? We now turn to a discussion of this aspect of the subject. First, we notice the 
following elementary 
NECESSARY CONDITION. Suppose F = (F,,...,F,) is a conservative vector field de- 
fined on an open set D C R", and suppose F is of class C'. Then 

OF} OF; 
Ons 7 Ox; 





for all i,j. 








Green’s theorem 15 


PROOF. We use the third characterization that F = Vf for some potential function f of 
class C!. Since F is assumed to be of class C!, we conclude that f is of class C?. Therefore, 
its mixed partial derivatives are equal, and we have immediately 


OF, . .d fof 

Ox; i Ox; (34) 
ee 
Ox; & 
_ OF; 
Ox; 











QED 


DEFINITION. Any vector field satisfying the condition we have just given is said to have 
zero curl. Also it is said to be irrotational. These terms will be explained in Chapter 13. In 
particular, we do not at the present time define the term “curl.” 

This necessary condition is certainly simply to apply. For instance, on R? the vector field 














(y?z + yz”, Qeyz+a2z7, xy? + xyz) 


is not conservative, since 


0 
py 2rye + £2") = Qry + 2xz, 
eZ 


0» 5 
—(xry + xyz) = 2ay+ x2 
By | y” + ryz) y 
are not equal. 
However, the necessary condition is not sufficient, even in case F is of class C!. Everyone’s 
favorite counterexample is the following. Use the polar coordinate “function” 6 on R?: 





r cos 8, 


x 


y= rang, 


Of course, @ is undefined at the origin and is also determined only up to additive integer 
multiples of 27. But the gradient of @ is completely well defined on R? — {0}. Using 6 = 
arctan(y/x) for example, we have 


Vi= : 
(ao a) 
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And that gives us our counterexample. Namely, the vector field 


is of class C®° on R* — {0}, and 

















OF, OF, 


‘Oy Ox’ 


but F’ is not conservative. For instance F' fails to satisfy the criterion 2, as we can detect if 
we integrate around any circle centered at 0: x = acost, y = asint for 0 < t < 27 gives 


QT : 
sint cost 
| Fed = | (-=. =) e (—asint, acos t)dt 
circle 0 a a 


20 
= | (sin? t + cos? t)dt 
0 


= 27. 


(Of course, we might formally write this computation as 


i. Fedt = | VO edz 
circle circle 


final point 


=0 





initial point 


= 2G) 


We shall call this nonconservative vector field V6; the vector field is well defined through the 
“potential function @” is not. 


PROBLEM 12-6. Do the differentiations to check that for the vector field V@ we 
have 

Oy Ox 
And then explain why it is completely unnecessary to do so, that the equation actually 
holds “automatically.” 
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PROBLEM 12-7. Let 7 be a closed curve in R? which loops around the origin exactly 
once in the counterclockwise direction: 
Show that ¥ 


i —ydz + xdy _ 
y 


27. 
x+y? 


PROBLEM 12-8. Determine all possible values of the line integral 
/ —ydzx + xdy 
y 


x+y? 


for curves y which start at (—3,0) and end at (1,—1) and do not pass through the origin. 


PROBLEM 12-9. Show that the vector field on R? — {0} given by 
3 2 
y —xty 
F(a, y) = (4 ) ’ 


rt? 


where of course r? = x? + y”, has zero curl. 


PROBLEM 12-10. Show that the vector field of the preceding problem can be 
expressed in the form 
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PROBLEM 12-11. Calculate the line integral 
/ y°dx — xy*dy 
» (2 +y)? 7 
where 7 is the counterclockwise ellipse “> — ve = 1. 


F. Sufficiency 





We now want to consider a C' vector field on an open subset D of R” and investigate 
whether it is conservative. Of course, we necessarily assume it has zero curl. We then turn to 
the question of deciding whether it is conservative. We’ll give two types of answers. 

CASE 1. n= 2, and D is simply connected. 


The hypothesis indicates that D has no “holes.” 
not simply connected 





simply connected 
Then it follows that zero curl = > conservative. We can see this from Green’s theorem. We 
verify criterion 2 of Section E. Consider a loop in D, and assume it is of a simple enough 
nature that it can be realized as the boundary of a nice region R: 





R 


Then Green = 


bdR bdR 


“fe 
= / | Odedy 


= 0. 
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This proves that criterion 2 holds for F’, at least for a pretty broad variety of loops. But 
certainly enough loops to ensure that the proof of 2=—>1= 3 holds. 


CASE 2. Any n, and D is all of R” or is a rectangle or .... 
We can go directly to determining a potential function by the simple device of integration. 
Before giving the proof, here is an example of the procedure. Let us take n = 3 and 





F = (2xy — 3yz, x? —3xz, 62? — 3zy). 


Without checking its curl let us just see whether there may be a potential function f(z, y, z). 
This would have to satisfy 

f. c= 2ry = 3Yyz ’ 

ii x? — 3rz 


fe 


627 — 3ry. 


Naively integrate the first equation: 
f =a°y — 3xyz + gly, 2), 
where g is a “constant” of integration. Plug this into the second equation to get 
2 ue 
x — 322+ Gy =X" — 3x2. 
Thus 
Gy = 9. 
Naively integrate this last equation to get 
g = h(z), 
where h is the ‘constant” of integration. Thus we now have 
f=2x°y — 8ryz + h(z). 
Plug this into the third equation to get 
—3xy + h'(z) = 62? — 3ay. 


Thus 
A@y= be. 


Integrate to get 
h(z) = 2z* + const. 


20 Chapter 12 


Thus we have succeeded in finding 


f =az°y — Bayz + 22°. 





MORAL. This naive procedure always works in case D is a rectangle or all of R”. If F is 
conservative, it produces a potential function. If F’ is not conservative, somewhere along the 
way an impasse will be reached and we'll learn that F is not conservative. 

Moreover, the general case works in exactly the same way as the example! We can proceed 
by induction on n. 

As an example, let us use the procedure to test whether the vector field 


F = (Qry — 3yz, 2? + 3xz, 62” — 3ry) 
is conservative. We suppose it is, so F' = Vf. Then the equation for Of /Ox gives as before 
f =2°y — 3ayz + gly, 2). 


Then the equation for Of /Oy gives 


O 
a? — 3024+ 4 =a? + 32z, 
Oy 
so that 
Og 
— = —-622Z. 
Oy 


This equation is impossible, as g = g(y,z) is independent of x, but the right side —6zxz is 
dependent on x. Thus F' is not conservative. 

Here is a description of the inductive procedure. We suppose F' = (F\,...,F,) is a Ct 
vector field defined on R", and we suppose F' satisfies the “zero curl condition.” That is, 























OF, OF; — 
= for 1<4,7 <n. 
Ox; Ox; ain 
Assume that we have verified the validity of our procedure for dimension n—1. We then apply 
it to the vector field (F,,...,/,-1) on R"!, where z,, is just regarded as a parameter. The 
integration procedure therefore produces a function fo(21,.-.,2%n—1,%n) which satisfies 
Ofo 


=f; for L<i<n-l. 





Ox; 


21 
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We anticipate that one more integration will give the desired potential function f. So we hope 
that 
f(x) = fo(x) + g(@n) 


will do. The required equation is Of /Ox, = F,, or in other words 


Ofo 
‘(t,) = F, - ——. 
g' (an) Dz, 
We need to know that the right side of this equation is independent of x,...,2%,_1. So we 


check its partial derivative with respect to x; for 7 < n: 


O(Fn— 92) OF, fo 














Ox; Ox;  Ox;OTp, 
eG O* fo 
Ox; 7 OF O08; 
Ox; 7 OX 


by the inductive assumption. Now the right side is indeed 0, thanks to the curl assumption. 
So one final integration gives us g(z,,), and, thus, f(z). 


G. A slick proof of sufficiency 


We suppose that F is a C' vector field on R”, defined on an open set B C R” which is 
star shaped with respect to one of its points 7g. This means that for every x € B, the line 
segment from 29 to x is contained in B. 

















star shaped convex not star shaped 


Any convex open set is star shaped with respect to any of its points. 


In what follows we assume (WLOG) that xo = 0. 
If F' is conservative with potential f, then from the equation Vf = F' we can derive a 
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formula for f. Namely, the FTC gives 
0. 
fe) =F) + | Sf(te)at 
0 
1 
=7(0) +f (Vf) (tx)  xdt 
0 
1 
~ #(0) + | F (tx) « xdt. 
0 
Now turn this around. Suppose F' satisfies the zero curl condition and define 
1 
ia= | F (ta) e xdt. 
0 


Notice that f is well defined on B, whether or not F' satisfies the extra condition; but the 
condition is exactly what gives us Vf = F’. We check this by performing a differentiation of 
f “under the integral sign” as follows: 


1 
a = a (| F(tx)e vdt) 
O 
I Ui 
O 


F) 
= ll = (F(tx)) ex + F(tz) a] dt, 











(F'(tax) e x)dt 





the last equality coming from the product rule. Now the chain rule gives 








a (F(te)) ex = (= ) ies 


“. OF; 
= ty Dn, (na: 


j=l 





nm 
=i 
= 


OF; 
Ox; 





Reg ee 


the last equality coming from the curl condition. The chain rule now gives 


0 0 
Fp (Plte)) #2 = t (Alta), 
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so the formula for ae becomes 








si = ; fe atta) + Fite) dt 
-{ © (tR(tx)] at 
=F, (tz) = 
= F;(x). 


Thus Vf = F. 


